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Abstract 

O 

Cosmological tensor perturbations equations are derived for Hamiltonian cosmol- 
^Jy ogy based on Ashtekar's formulation of general relativity, including typical quantum 

gravity effects in the Hamiltonian constraint as they are expected from loop quan- 
tum gravity. This translates to corrections of the dispersion relation for gravitational 
waves. The main application here is the preservation of causality which is shown to 
! be realized due to the absence of anomalies in the effective constraint algebra used. 

CO 

: 1 Introduction 

o . 

In cosmology, the study of gravitational waves created through physical processes in the 
early universe provides a unique window to initial stages of the universe. Significant efforts 
are being made to detect possible signatures of tensor mode perturbations of space-time 
geometry through measurements of the polarization in the Cosmic Microwave Background 
(CMB). Since quantum gravity effects could play a significant role in the very early universe, 
it is of interest to study possible quantum gravity effects on gravitational wave propagation 
during these periods. In the last few years, applications of the method used in loop 
quantum gravity (LQG) pQElE], a candidate quantum theory of gravity, to early universe 
cosmology have led to significant progress. In particular, the quantization of homogeneous 
cosmological models known as loop quantum cosmology (LQC) [I] has led to a resolution 
of the big bang singularity [HI El El El E] , an d techniques have become available to include 
inhomogeneous perturbations [TOT ITT]. 

In this paper we study typical quantum gravity effects for tensor modes that are ex- 
pected from loop quantum gravity^ In particular, we consider the effects on gravitational 
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1 Possible cosmological implications for primordial gravitational waves with input from loop quantum 
cosmology have recently been discussed in [12] . However, only corrections to the background dynamics were 
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wave dynamics expected from corrections to classically divergent inverse powers of metric 
components and from the use of holonomies in the quantum theory instead of connection 
components. To study the dynamics we compute gravitational wave equations together 
with their dispersion relations. 

In section [3j we present a derivation of tensor mode equations in Ashtekar variables. 
The calculations are purely canonical and split off the tensor mode in the metric from the 
outset. This mimics the usual covariant derivations as far as possible in a way accessible to 
canonical quantizations. Other canonical derivations exist [13] which due to their explicit 
use of Dirac observables appear more difficult to use in quantizations. In the following 
section, we consider the effects of quantum corrections to the inverse volume in Hamil- 
tonian and compute the correspondingly corrected tensor mode equation. After that, we 
consider a second quantum effect due to the use of holonomies in a loop quantization. From 
the corrected wave equations one can easily derive the corresponding dispersion relations. 
Quantum gravity corrections are sensitive to the underlying discreteness of a quantum 
state, which in general changes as the universe expands. Thus, also propagation speeds de- 
rived from the corrected dispersion relations are functions of time, providing in particular 
a varying speed of light scenario. 

Both corrections are typical of loop quantum gravity and thus test its basic features. 
The dispersion relations in particular allow one to investigate possible violations of causal- 
ity which would arise if the propagation velocity of gravitational waves would turn out to 
be larger than the speed of light. In fact, we will see that gravitational waves travel faster 
than the classical speed of light, but not faster than the physical speed of light which is 
also subject to quantum corrections from an underlying discrete geometry. Quantum cor- 
rections to the gravitational and electromagnetic dynamics are related by the requirement 
of anomaly freedom, which can directly be implemented at the effective level and implies 
that physical causality is preserved. 

2 Canonical formulation 

We consider linear tensor mode perturbations around spatially flat Friedman-Robertson- 
Walker (FRW) spacetimesJl The general form of a perturbed metric around the isotropic 
background FRW background containing only the tensor mode is 

g Q0 = -N 2 + q ab N a N b = -a 2 ; g 0a = q ab N b = ; 9ab = q ab = a 2 [5 ab + h ab ] , (1) 

considered while perturbation equations for the tensor mode were otherwise left unchanged. Moreover, 
the main analysis there focuses on non-perturbative quantum effects in the background dynamics which 
makes a rigorous inclusion of perturbative inhomogeneities around the background difficult; see e.g. the 
discussion in [10]. In the present paper, by contrast, we provide a consistent perturbative setting in which 
quantum corrections of the inhomogeneities themselves arc included in the equations. As we will see, this 
by itself provides important effects which are not mimicked by corrections to the background dynamics. 

2 The procedure follows that used for scalar [TT] and vector modes [2] but is simpler at several places 
in the derivation of their equations of motion as well as for gauge issues. 
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where a(t) is the scale factor of the FRW spacetime. This notation is adapted to a canonical 
formulation, where the spacetime metric is decomposed in terms of the spatial metric 
q a b, the lapse function iV and the shift vector N a . Here we use the convention that Greek 
letters denote space-time indices whereas small Latin letters denote spatial indices. The 
symmetric metric perturbation field h a b is transverse and traceless, i.e. it satisfies d a h a b = 
and 5 ab h ab = 0. This removes any vectorial or scalar contributions from gradient terms 
<9( a t>ft) or d a dbV or from the trace uS a b which rather contribute to the vector and scalar 
modes. Also the lapse N and shift N a , being scalar and vectorial, respectively, do not 
contribute to tensor perturbations. Thus, in a canonical formulation tensor perturbations 
are generated through perturbations of the spatial metric q a b alone. 

2.1 Background 

In Ashtekar's formulation of general relativity [H2, [16], the spatial metric as a canonical 
field is replaced by the densitized triad Ef, defined as 

£?:=|det(ej)|e?. (2) 

Here, el as a matrix is the inverse of the co-triad e l a whose relation to the spatial metric 
is q a b = e*e|. The canonically conjugate variable to the densitized triad is the Ashtekar 
connection A l a := T l a + jK^, where K l a is the extrinsic curvature and 7 is the so-called 
Barbero-Immirzi parameter [TT1 [T6] . The spin connection T l a is defined such that it leaves 
the triad covariantly constant and has the explicit form 

K = -e^e){d [a et ] + \ele l a d [c e l b] ). (3) 

As we perturb basic variables around a spatially flat FRW background, our background 
variables denoted by a bar are 

Ei = pS? ; ^ = ; K = W a ; N = VP ! N a = , (4) 

where p = a 2 and the spatial metric is q a b = a 2 #afcH The choice of iV = a leads to conformal 
time which is used in what follows. 

2.2 Perturbed canonical variables 

The perturbed densitized triad Ef and Ashtekar connection A % a around a spatially flat 
background are given by 

Et = p6? + 5Et ; A\ = P a + jK^ = + (6V a + ^K) , (5) 

where p and 7^ are the background densitized triad and Ashtekar connection, using the fact 
that f = for a spatially flat isotropic model. The general form of a co-triad corresponding 

3 Compared to [18j we drop an additional tilde on p to keep the notation simple. 
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to a spatial metric as in ([I]) is 



el 



51 + \K 



(6) 



where h\ := 5 lb h a b. The densitized triad ([2]) then has the perturbation 



5E? = - l -ph1 , (7) 

where we have used the fact that tensor mode perturbations are traceless, i.e. S l a 5Ef = 0. 
For a general perturbed densitized triad (jSJ) the linearized spin connection (jHJ) becomes 

SVi = -JHa&SE* ■ (8) 
p J 

As perturbations of lapse N and shift N a do not contribute to the tensor mode, we can 
set 5N = and SN a = when studying tensor mode dynamics. 

As described in more detail for scalar and vector modes in [TTJ, [F4], the symplectic 
structure splits into one for the background variables and one for perturbations, 

fk,P} = ^ , {5K(x),SE b 1 (y)} = 8nG5 3 (x,y)5 b a S l J . (9) 

Here, G is the gravitational constant and Vq is a fiducial volume introduced to arrive at a 
finite symplectic structure for the background variables by integrating the action only over 
a finite cell rather than all of IR 3 . Since this background is homogeneous, no information 
is lost by the restriction to a cell. However, a fiducial quantity enters the formalism which 
must disappear from final physical results. 

This provides separate canonical structures for the background and perturbations, but 
these variables will be coupled dynamically. In particular, the homogeneous background 
dynamics would receive back-reaction effects at quadratic or higher order. 



3 Classical dynamics 

In canonical quantum gravity, dynamics is determined by a Hamiltonian (constraint) oper- 
ator rather than a path integral. This implies, that one obtains relevant quantum correc- 
tions at the level of an effective Hamiltonian as opposed to an effective action in a covariant 
quantization. To study the effects of quantum corrections to the classical equations of mo- 
tion one thus needs to derive these equations starting from an effective Hamiltonian. Here, 
we are interested in studying the quantum correction expected from loop quantum gravity 
which is based on Ashtekar variables in the classical formulation. As a preparation for an 
analysis of effective tensor mode Hamiltonian we thus derive in this section the classical 
gravitational wave equation in canonical gravity using Ashtekar variables. 
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In a canonical triad formulation of general relativity there are three types of constraints: 
the Gauss constraint which generates local rotations of the triad, the diffeomorphism con- 
straint which generates spatial diffeomorphisms and the Hamiltonian constraint which 
completes the space-time diffeomorphisms and is thus relevant for the dynamics. For lin- 
ear perturbations including only the tensor mode, the corresponding Gauss constraint is 
trivially satisfied as the perturbation field h % a = 5 lb h a i, is symmetric. In fact, the triad 
perturbation ([6]) is symmetric, while su(2)-gauge transformations of the Gauss constraint 
could only generate antisymmetric contributions owing to the antisymmetry of the sub- 
structure constants. Also the diffeomorphism constraint is identically satisfied as N a = 
for the tensor mode as discussed before. Thus, solutions for tensor mode perturbations are 
completely governed by the Hamiltonian constraint. 



3.1 Hamiltonian constraint 



The Hamiltonian constraint generates 'time evolution' of the spatial manifold in terms of 
a time coordinate. Its general expression is 



H G [N] 



i6ttG y E 



d xN 



v/|det£| 



- 2(1 + 1 2 )K{K^ 



(10) 



Using the expression (jSJ) of the perturbed basic variables and the curvature F^ b = d a A b 



dbA l a + eijkA{A b , one can simplify and expand (TIDl) for linearized tensor modes. Up to 
quadratic terms in perturbations we have 



H G [N] 



16nG 



d xN 



-Qk 2 ^p- 



2k 1 



As expected, 7 dependent terms drop out of the Hamiltonian constraint when one uses the 
spin connection and the fact that densitized triad and extrinsic curvature are symmetric 
for the tensor mode. 



3.2 Linearized equations 

In the standard covariant formulation linearized equations for metric perturbations are 
derived by considering the variation of the action with respect to the perturbed metric. In 
a canonical formulation, the linearized equations are derived using Hamilton's equations 
of motion. For the perturbed densitized triad, 

SEf = {6E%, Hq[N] + H matter [N]} (12) 

leads to the expression of extrinsic curvature. Here H matteT [N] denotes the matter Hamilto- 
nian which together with the gravitational contribution (TTTf) forms the total Hamiltonian. 
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Also the matter Hamiltonian depends on the lapse function through the determinant of 
the space-time metric. The choice of the background lapse function then determines the 
time coordinate which the dot refers to, which from now on will be N = a for conformal 
time. 

Using the expression (J7J) of the perturbed densitized triad and thus SE® = — ^(ph? + 
phf), the equation of motion ffT2l) then leads to the expression 



5Kj 



K + kh\ 



(13) 



for the linearized extrinsic curvature, where we used the background extrinsic curvature 
k = p/2p which follows in a similar way from the zero order Hamiltonian constraint. 
The second Hamilton equation of motion 



SK a = {SK l a , H G [N] + H mMer [N}} 



(14) 



describes the evolution of perturbed extrinsic curvature. Using (fT3l) . one can derive the 
second order equation of motion for gravitational tensor mode perturbations: 



hi + 2kK - v 2 K 



8ttGTP 



where 



it. 



1 dH, 



matter 



SEjfilfil 



+ 



5H n 



(15) 



(16) 



W dp V P J ' 5 ( 5E i) 
As usually, in the absence of source terms (115p has propagating wave solutions which are the 
usual gravitational waves in the given cosmological background. Cosmological expansion 
leads to a friction term which is proportional to k and thus the Hubble parameter. 

The quantity describes the linear transverse and traceless source terms that can 
be related to the transverse and traceless part of the perturbed stress-energy tensor as 
IT; = p6T® l a . For comparison, we now demonstrate the explicit relation between and 
the stress-energy tensor which is defined as 



2 5S- 



matter 



-9 $g 



(17) 



for a given matter action Scatter- Including only tensor perturbations, the inverse spatial 
metric can be written as q ab = g ab since the shift vector iV a does not contribute to tensor 
perturbations, i.e. iV a = 0. Thus, space-space components of the stress-energy tensor are 
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2 6H, 



matter 



-ab 



[N] 



-9 Sg 



ab 



ab 



where q is the determinant of the spatial metric 



The inverse spatial metric q is 



related to the densitized triad by qq 



ab 



E?E!?, for which we use the perturbed form 



Ef = El + bE^ with Ef = p5f. In perturbation theory both Ef and 8Ef are treated as 
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independent degrees of freedom and one can express the stress-energy tensor (fl8l) up to 
linear order in perturbations as 



T, 



ab 



-^matter 




6H, 



matter 



5(5E 



(19) 



With the relation between the inverse spatial metric and the densitized triad, one can show 
that 



.dE c - 
5 3 — 3 - 

"dq ab 



V 



SE C 



Sab + 



2p 



dq 



ab 



E c 



V 
~2 



[5^5 j)b -5]6 ab ] (20) 



using the fact that tensor perturbations are symmetric and trace-less i.e. S 3 c 5Ej = 0. While 
in the first equation we have kept terms up to first order in perturbations, in the second 
equation we have kept only the zeroth order terms as the term (5H m /5(5Ej)) itself is at 
least of first order in perturbations. We then compute the perturbed stress-energy tensor 



5T a :=T a -T a 



Sl(q cb T c 



ab) 



rpl 



1 dH matter ^E^5 l c 6i) 5H matteT 

+ 



V 



5{5Ef) 



-= IT 

N^p a ' 

(21) 

where we have used the requirement that for tensor perturbation, perturbed stress-energy 
tensor is trace-free, i.e. 5j(5H mautteT /5(5Ej)) = 0. The background stress-energy tensor T % a 
is given by 



61 



NV p 3 / 2 



2pdH, 



matter 



dp 



[22] 



This expression explicitly shows the relation between spatial components of the background 
stress-energy tensor and background pressure. 



4 Quantum dynamics 

In the previous section, we have seen how the tensor mode equation is derived from canon- 
ical classical cosmology. We will now include two basic types of quantum corrections that 
are expected from the Hamiltonian of loop quantum gravity. These corrections arise for 
inverse powers of the densitized triad, which when quantized becomes an operator with 
zero in the discrete part of its spectrum thus lacking a direct inverse [19] , and from the fact 
that a loop quantization is based on holonomies, i.e. exponentials of the connection rather 
than direct connection components. There is an additional source of corrections due to 
back-reaction effects of quantum fluctuations on expectation values of the basic variables 
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[20| [2T] . This is more complicated to derive and not included in the present analysis. We 
need to consider these corrections only in the Hamiltonian constraint because the full dif- 
feomorphism constraint does not receive quantum corrections. It thus remains trivial for 
the tensor mode dynamics as in the classical case. 



4.1 Inverse volume corrections 



In loop quantum gravity, the factor EjE%/ y^det E\, which appears in the Hamiltonian 
constraint ( fTOl) and contains inverse powers of the densitized triad, cannot be quantized 
directly but only after it is re-expressed as a Poisson bracket not involving an inverse [19] . In 
homogeneous models, explicit calculations show that eigenvalues of the resulting operator 
approximate the classical expression for large values of densitized triad components, but do 
provide quantum corrections which become larger for small components [22l 1231 [24] . One 
can include these corrections as one of the new terms in effective expressions by introducing 
a factor a whose generic form in the large volume regime is 



a(p) = l + c(^l , (23) 

where n and c are positive numbers. Anticipating similar quantum corrections even for the 
inhomogeneous case, the effects of such a correction have already been studied for scalar 
and vector mode perturbations [TTj [T4") . Here, we provide an analysis for tensor mode 
perturbations, starting with a corrected Hamiltonian constraint 

Ht n [N] = J^l d 3 xNa(p, 5E?) [-6fc 2 yf - ^(SE^E^) + VP(5Ki5K k d 5^) 

2k 1 
- ^E]SKi) + ^(6 cd 5^d e Ep f E d k )^ .(24) 

(We indicate quantum corrected expressions by a superscript "phen" to indicate that such 
terms are introduced for a phenomenological analysis while a systematic effective analysis 
is still outstanding.) This is to be used in a perturbative inhomogeneous context and is 
thus not set in a purely minisuperspace model. In this case, a(p, SEf) also depends on 
triad perturbations and is in general more complicated to compute from an underlying 
Hamiltonian operator than in homogeneous models. Moreover, since the function a comes 
from the quantized inverse densitized triad where the tensorial term EjE^/ a/| det E\ is 
quantized as a whole, it could be tensorial in nature. However, later we will see that its 
leading effect on perturbation dynamics comes from the background corrections a(p, SE® = 
0) = a. 

The only background variable determining the geometry is p, as a function of which the 
corrections are expressed. The appearance of such a scale factor dependent function in dy- 
namical equations has occasionally led to concerns that quantum gravity might break the 
scale invariance of flat isotropic models, or even introduce gauge artefacts. Alternatively, 
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one can absorb the rescaling freedom in a redefinition of the fiducial volume Vo encoun- 
tered earlier, but then the dynamical equations as well as their solutions seem to depend 
on this fiducial volume. None of these problems occurs in genuine inhomogeneous models. 
The dependence of a correction function a in an inhomogeneous Hamiltonian constraint is 
through elementary area variables whose values are determined by an underlying inhomo- 
geneous state. (Areas, or more precisely fluxes, are elementary because they are directly 
related to the densitized triad as a canonical variable.) These elementary areas build up 
the quantum geometry of space in a discrete manner and their sizes determine the degree 
of discreteness involved. The scale of corrections, too, is determined by the underlying 
state and thus depends on the size of discreteness. 

More specifically, correction functions only seem to depend directly on the scale factor a 
because other parameters, most importantly the number of lattice sites Af per volume in the 
underlying state, have been suppressed (see also [25]). This parameter rescales in the same 
way as the scale factor such that the whole expression is scaling invariant. Elementary areas 
are the primary object appearing in corrections and they are, on average, of the geometrical 
size F = a 2 £l where £o is the average coordinate length of lattice links. This quantity is 
certainly scaling-independent. Moreover, £ is related to Af and thus depends on the precise 
quantum state and has to be determined from the underlying theory. The parameter Af, 
however, also depends on the chosen volume Vo in which one counts the number of lattice 
sites: Af = Vq/1$. One can identify Vo with the fiducial volume introduced earlier. Then, 
an alternative worry has been voiced, namely that a scaling invariant quantum correction 
would depend explicitly on the fiducial volume. Also this is not true: One simply rewrites 
the quantity F as before in a different way, F = a 2 V^ 3 /Af 2 / 3 . Numerator and denominator 
are now scaling independent but Vo-dependent. Nevertheless, the total quantity F which 
appears in quantum corrections from the inverse volume is Vo-independent. Such quantum 
corrections are thus consistent and do not depend on any gauge or other choices. 

4.1.1 Linearized equation 

Extrinsic curvature is derived using Hamilton's equation of motion (jl2"]) . With quantum 
corrections in the Hamiltonian, also extrinsic curvature should receive quantum corrections. 
In our case of a Hamiltonian (T24"j) . this leads to 



Here one can see that the leading correction due to the background correction function is a 
as inhomogeneous contributions to a will contribute only higher order terms. The second 
Hamilton's equation together with the just derived expression of extrinsic curvature (|25|) 
then provides a second order equation 



(25) 



1 \-h l a + 2k(l- —1 hi - aV 2 ht\ + A\ = 8irGlli 



a \ a J 



(26) 
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for the dynamics of tensor mode perturbations, where the prime denotes a derivative by p 
and 



A\ = 3Nk 2 y/p 



_d{8Ef) 3pdp 

Inverse densitized triad corrections lead to several significant changes in the wave equa- 
tion (1261) compared to its classical counterpart (|15p . First, there are corrections in the 
coefficient of h l a and the coefficient of the Laplacian term V 2 h l a . Secondly, there are ad- 
ditional contributions to the friction term and, thirdly, an entirely new term A l a . In the 
context of vector mode dynamics [13] , the same term A l a appears in the equation of mo- 
tion but it also presents an anomaly term in the constraint algebra between the perturbed 
Hamiltonian and diffeomorphism constraints. Requiring an anomaly-free constraint alge- 
bra in the presence of quantum corrections then implies that A l a must vanish and leads to 
restrictions on the possible functional form of the quantum correction function a(p, 5Ef). 
While there are no such anomalies in the constraint algebra for tensor modes as the diffeo- 
morphism constraint is trivial here, the same quantum correction function a(p, 5E®) as for 
the vector mode must occur since there is only one Hamiltonian constraint which is just 
split into different mode contributions to simplify the analysis. Thus, we must set A\ to 
zero, which we will do in the subsequent analysis. 



4.2 Holonomy corrections 

A loop quantization represents holonomies as basic operators on a Hilbert space rather 
than connection components. Moreover, it is impossible to derive operators for connec- 
tion components from holonomies and thus any quantized expression depending on the 
connection must do so through holonomies. This is especially true for the Hamiltonian 
constraint, which thus receives quantum corrections from higher powers of the connection. 
Holonomies are non-linear as well as (spatially) non-local in connection components. Thus, 
they provide higher order and higher spatial derivative terms. Higher time derivatives, as 
they would also be provided by higher curvature terms, do not arise in this way but rather 
through the coupling of fluctuations and higher moments of a quantum state to the expec- 
tation values [2U1 EI]. Here we focus on corrections from holonomies as a typical effect of a 
loop quantization, while the more complicated quantum back-reaction effects are genuine 
and occur for any interacting quantum theory. 

We start by recalling the situation for a homogeneous and isotropic model with a 
massless free scalar field. This model allows one to compute an exact effective Hamiltonian 



H$[N) 



16nG 



f28) 



where higher order terms of extrinsic curvature (which is proportional to the Ashtekar 
connection in a spatially flat model) are explicit in the sine. Again, Vq is the volume 
of the fiducial cell introduced to avoid the integration over spatial infinity in (|T0l) for a 
homogeneous background. Moreover, p, is a new parameter related to the action of the 
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fundamental Hamiltonian on a lattice state. It can be understood as the coordinate size 
of a loop whose holonomy is used to quantized the Ashtekar curvature components F l ab . In 
the limit p, — > 0, the effective Hamiltonian reduces to the standard classical Hamiltonian. 
In general, p can even depend on the triad component p to reflect refinements of the 
discrete state during dynamics [25] . While the precise behavior is difficult to compute, 
general considerations restrict the dependence to p(p) = where < n < — 1/2. Only 
the limiting cases n = [271 [18] and n = —1/2 [28] have so far been discussed in the 
literature. 

Variation of the Hamiltonian constraint with respect to the background lapse function 
TV leads to the effective Friedmann equation 



— 2 

1 / sin/17/A 8tvG 



V \ PI 

where p is the energy density defined as 



P , (29) 



5H, 



P 



matter 



(V p 3 / 2 ) 5N 



In this form, the effective equation is precise only for the energy density of a free scalar, 
-^matter = \NVop~ 3 l 2 p 2 i i with momentum p^. If a matter potential or anisotropics and 
inhomogeneities are added, additional corrections arise [21] from quantum back-reaction. 

While classical cosmological dynamics is in general singular, the effective dynamics 
is non-singular. The singularity avoidance is achieved by exhibiting a bounce at small 
volume when the energy density reaches a critical value [8]. This can be seen by writing 
the effective Friedmann equation (1291) as (smp'jk) 2 = p/p c where 

3 , , 

P c = a n-2 2- ■ ( 30 ) 

The boundedness of the sine then implies a minimum p and thus a minimum non-zero 
volume, the bounce scale. For the case p = y^KJp for instance, p c is a constant. The 
critical energy density p c then signifies the maximum energy density that is reached at the 
bounce point. This can be seen explicitly from the Hamilton's equations of motion which 
are 

p = 2 NVp( S ^) (31) 



2fl J 



and 

op 



Vp 



sin p^k^ 
PI 



]7TG dH ma tter ^2 -j 



W dp 



Thus, for an isotropic model sourced by a massless, free scalar field the effective Hamil- 
tonian can be obtained by simply replacing the background Ashtekar connection 'jk by 
p~ x sin p'yk. This is no longer true for other models, especially when inhomogeneities are 
included. But to study the effects on inhomogeneous perturbations, one can substitute 
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the appearance of k in the classical Hamiltonian by a general form smmM7fc where m is a 
number. There may be additional corrections, but qualitative effects can already be read 
off from such a replacement. With this prescription, the Hamiltonian constraint becomes 



Hi 



phen 



[N] 



16nG 



d 3 xN 



6 ^/sin^fc 
V ^7 



2 p3/2 



sin pirjh 

Pn 



{8E c JEi8%) 



+ ^p{5Ki5K k d 5 c k 5" 



+ 



3/2 



J7 



/A- 



8 ef d e Ep f E%) 



, (33) 



In writing the explicit coefficients we have required that the Hamiltonian has a 'homoge- 
neous' limit in agreement with what has been used in isotropic models ( 128]) . This fixes 
the parameter m to equal one in the first two terms. The parameter for the last term as 
chosen here is the one which leads to an anomaly-free constraint algebra in the context of 
vector modes [H]. 

One should keep in mind that, although we write explicit sines in this expression and 
thus arbitrarily high powers of curvature components, this is to be understood only as a 
short form to write the leading order corrections. This is more compact than writing the 
leading terms of a Taylor expansion of the sines. The expressions are, however, reliable 
only when the argument of the sines is small, which excludes the bounce phase itself. 
Moreover, higher orders are supplemented by further, yet to be computed higher curvature 
quantum corrections. (Such sine corrections can be used throughout the bounce phase 
only for exactly isotropic models sourced by a free, massless scalar [2"B"l [SO] .) 



4.2.1 Linearized equation 

The expression for extrinsic curvature is again derived using one of Hamilton's equations 
of motion and thus receives quantum corrections also from the use of holonomies in loop 
quantum gravity, 



SK n 



hi + 



f sin 2^jk\ 



(34) 



V 2>y ; 

Along with Hamilton's equation for the perturbed extrinsic curvature, this equation (I34|) 
then yields the quantum corrected second order equation for tensor perturbations 



V rl J 



(35) 



This equation ( 1351) describes propagating degrees of freedom which are the usual gravita- 
tional waves subject to quantum corrections. Unlike for inverse densitized triad corrections, 
the coefficients of h l a and V 2 h l a take the classical form. On the other hand, the friction 
term does receive corrections. As a new feature, there is an additional term proportional 
to field perturbations h\ with coefficient 



T, 



Q 



p dfi 
ft dp 



-2 2 

/x 7 



sin fijk 



(36) 
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For any /i oc |p| n with n < 0, Tq is positive definite. Finally, the source terms from the 
matter Hamiltonian take the form 



n, 



Qa 



1 OHr, 



3V dp 



P 



cos 2/ryfc + 



5H n 



5(5E a i 



(37) 



as the transverse and traceless part of the stress-energy tensor that sources gravitational 
waves. The additional cosine can be understood from the fact that the background geom- 
etry receives quantum corrections and is used to define the trace-free part of stress-energy. 
The source Hq 1 vanishes when there is no matter field, and it reduces to the classical 
transverse and traceless part of the stress-energy tensor 11^ in the limit p, — > 0. 



5 Dispersion relation 

To study wave propagation it is often convenient to compute the relevant dispersion relation 
from the corresponding wave equation, presenting a relation between the frequency and 
the wave vector. In this section, we use dispersion relations for the quantum corrected 
gravitational wave equations to study some of their basic properties. 

Starting with the classical dispersion relation to be able to contrast it with the corrected 
versions later on, we consider the source-free tensor mode perturbation equation by making 
a plane wave ansatz h l a oc h l a exp(iut — ik.x). Here, the frequency u corresponds to proper 
time t where the lapse function N, in contrast to the previous section, is equal to unity. 
The classical tensor mode equation (fl5|) then simply implies 




Here we have ignored the friction term in the equation of motion (fT5|) since we are mainly 
interested in local propagation not involving cosmic scales. The dispersion relation (138]) is, 
of course, the standard classical dispersion relation between the frequency u and the proper 
wave number k/a. We further note that the corresponding group velocity of gravitational 

waves 

duj 

Ugw := dW^y 

is equal to 1 (in natural units). 



5.1 Inverse volume corrections 

Repeating the calculations of the classical case but using quantum corrections in the wave 
equation we obtain the corrected dispersion relations. In particular, the tensor mode 
equation (1261) in the presence of inverse volume corrections leads to 
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Dispersion relation with inverse volume corrections 



Classical 
Corrected (a-,) 
Corrected joy 
Corrected (a 3 ) 



(a-f > a 2 > 0L3) 




Proper wave number (k/a) - 



Figure 1: Dispersion relation for gravitational waves in the presence of inverse volume 
corrections. 



as illustrated in Fig. [IJ 

As one can see, the quantum correction function multiplies the wave number k, thus 
affecting the mode on all scales. Moreover, given that a > 1, the corrected group velocity 
due to inverse volume corrections is greater than unity. This may appear as a violation 
of causality since gravitational waves would travel faster than with the speed of light. 
However, this refers to the classical speed of light, while a physical statement requires us 
to compare the velocity to the physical speed of light. This differs from the classical one 
because also the Maxwell Hamiltonian receives inverse volume corrections in loop quantum 
gravity [31]. In the regime of linear inhomogeneities such corrections have been computed 
in [32], and a derivation of the quantum corrected group velocity of electromagnetic waves, 
which we present in Sec. El shows that it is not smaller than that of gravitational waves. 
Thus, there are no violations of causality. 

5.2 Holonomy corrections 

We finally consider corrections to the dispersion relation of gravitational waves due to the 
appearance of holonomies. Again ignoring the friction term, a plane wave ansatz in the 
wave equation fl35l) leads to 

l)\ mi (40) 
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Dispersion relation with holonomy corrections 



Classical - 
Corrected (m .) - 
Corrected (m „) - 
Corrected (m ™J 

(m 1 g > m 2 g > m 3 g ) 




Proper wave number (k/a) — > 



Figure 2: Dispersion relation in the presence of holonomy corrections. The classical dis- 
persion relation is approached as the background energy density decreases. 



where 

One may note that holonomy corrections effectively contribute a new additive term m 2 
in the dispersion relation (j4TJl) compared to the classical dispersion relation fl38l) . With 
this quantum correction, the gravitational wave has acquired an 'effective mass'. This, in 
turn, implies that the different modes of the gravitational waves propagate with different 
group velocities which are less than unity. Also here, causality is thus respected because the 
curvature independent electromagnetic Hamiltonian does not receive holonomy corrections. 
The corrected dispersion relation (HUj) is shown in Fig. [2j 

Using expression (l4Tj) . one can estimate the value of the 'effective mass' of the graviton 
at the present epoch. Given the value of 7 ~ 0(1) , A ~ 0(l)£f>, p c ~ 0(1)M| and 
the energy density p ~ 10~ 120 Mp of the present universe one obtains the value m g ~ 
10 -120 Afp = 10 _92 eV. Here £p and Mp are Planck length and mass respectively. Current 
observational bounds on the graviton mass from solar system measurements is m g < 4.4 x 
10~ 22 eV and its accuracy could be lowered up to m g < 10~ 26 eV from future gravitational 
wave measurements [33j [3U [35l [36] . Thus, our estimated theoretical value of the 'effective 
graviton mass' is well below the observational bound at present. It is unlikely that such a 
value could be tested observationally in the near future. However, given that the 'effective 
mass' depends on the background energy density, such an effective mass could play a 
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significant role in early universe physical phenomena such as inflation. 

6 Causality 

To determine whether causality is respected by the quantum corrections, we have to com- 
pare the propagation speed of gravitational waves to the physical speed of light. Just as 
tensor perturbations of the metric receive quantum gravity corrections, the electromagnetic 
field also is corrected. Thus, the speed of its wave excitations may differ from the classical 
value in the same way in which the gravitational wave velocity differs from the classical 
one. For an analysis of causality the two corrected velocities have to be compared. 

The basic field of Maxwell's theory of electromagnetism is the vector potential A^. Its 
source-free dynamics in a general space-time background is governed by the action 



where the background space-time is specified by the Lorentzian space-time metric g^ v . We 
again use the convention where Greek letters denote space-time indices whereas small Latin 
letters denote spatial indices. 

6.1 Canonical formulation of the electromagnetic field 

In a canonical formulation, as before, the space-time metric is decomposed into the spatial 
metric q^ v = + n fl n u and normal components which provide the non-dynamical lapse 
function and shift vector. Also the electromagnetic fields are decomposed, with the electric 
field 7r a arising from the space-time components of the field strength, and the purely spatial 
components giving the magnetic field. The total Hamiltonian for the electromagnetic 
field corresponding to the action (|42p can be written as 




(42) 



7^em = Hem [N] + -Dem [N a ] + Gem 



(43) 



where Hem [N] denotes the electromagnetic contribution 




(44) 



to the Hamiltonian constraint. Similarly, there is a contribution 




(45) 



to the diffeomorphism constraint, and a U(l)-Gauss constraint 




(46) 



16 



As usually in canonical formulations, time and space components of the physical fields such 
as A Q and A a play different roles for the dynamics. Variation of Gem with respect to A , 
the time component of A^, whose conjugate momentum is absent in (1431) and which is thus 
a Lagrange multiplier, leads to the usual expression 

d a 7r a = (47) 

of the Gauss law. 

Hamilton's equations of motion for the canonical fields A a and 7r a take the form / = 
{/,Hem}, explicitly given by 

AttN 

A a = d a A + N c F ca + —n c q ca , (48) 

and 

tt° = d c (N c n a ) - d d (N a 7r d ) + ^-d c (N^F ef q ce q^) . (49) 
For further details of the canonical analysis we refer to [32] . 



6.2 Classical propagation 

As before, we analyze the propagation of linear electromagnetic waves on a spatially flat 
Friedmann-Robertson-Walker background. (Small perturbations of the electromagnetic 
wave will induce small perturbations for geometric variables as well. However, to linear 
order the perturbations are independent of each other and can thus be studied separately.) 
As before, the spatial metric q ab and shift vector N a then are 

q ab = a 2 5 ab ; N a = , (50) 

where a(t) is the scale factor. 

The momentum can be eliminated from (148]) by computing its divergence and using the 
Gauss constraint (j47p : 

d t (d a A a ) - V 2 A t = (51) 

where dt refers to the time derivative according to Hamilton's equations of motion, d c = 
5 ec d e and S ec d e d c = V 2 . To satisfy this equation we make the standard gauge choices A = 
and d a A a = 0. The equations (]4"8l) and (]4"9l) together then lead to the electromagnetic wave 
equation 

d t (^d t A a ) - V 2 A a = (52) 

With the choice of proper time i.e. N — 1, the wave equation fl52l) has the usual friction 
term due to the evolving cosmological background whereas with the choice of conformal 
time i.e. N = a, the friction term will not be explicitly present. 

From the equations of motion we can again compute the dispersion relation using the 
standard wave ansatz A a ~ A a expi(ut + k.x). Here, we will choose N = 1 so that the 
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frequency ui corresponds to proper time, but we ignore the friction term which is justified 
for small wave lengths compared to cosmological scales. The classical wave equation (|52|) 
then leads to the standard dispersion relation 

✓ =g)'. (53) 
Moreover, the group velocity of electromagnetic wave propagation is 

VEM = = ( } 

which is constant in a classical cosmological background. 

6.3 Propagation in the presence of quantum gravity corrections 

Quantum gravity corrections mainly affect the Hamiltonian constraint, which becomes 

E2UEHE2] 



aEM(qcd)^pTr a ir b qab + pEM(q c d)^-F ab F cd q ac q b 



(55) 



where a E M(q c d) and /?EM(<?cd) are the correction functions due to quantum gravity effects 
in inverse triad components. This provides equations of motion 

AttN 

A a = d a {t»A^) + N c F ca + -—a EM 7i c q ca , (56) 

V ^ 

and 

n a = d c (N c n a ) - d d {N a ir d ) + ^d c (Np EM ^F ef q ce q a f) . (57) 

Using the same gauge fixing A t = and d a A a = 0, which is possible since the Gauss 
constraint does not receive quantum corrections, one obtains the corrected wave equation 

a a „ \ ( N@EM \ ^2 



d t -r^z — d t A a - V 2 A a = , (5* 

Na EM J \ a 1 



where a EM ■= a EM \ qcd =^5 cd and (3 EM ■= l3 E yi\q cd=a H cd - This provides the dispersion relation 

uj 2 = a EM p EM (59) 

and group velocity 

V E M = j703 = V " EM ^ EM ' ( 60 ) 

As in the case of gravitational waves, we see that v E m > 1 is larger than the classical 
value, since both qem and /?em are always greater than one in perturbative regimes [TO] . 
Secondly, the group velocity is no longer constant but varies with time as the universe 
expands. Varying speed of light has been studied in literature mainly in the cosmological 
context [581 159"]. and motivated for instance from bi- metric gravity [40J or non-commutative 
geometry [UJ; see [12] for a review. 
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6.4 Relation between the speed of gravitational waves and the 
speed of light 

In the electromagnetic Hamiltonian fl55|) we have seen two quantum correction functions 
qem and /?em- For homogeneous situations, they have the generic feature of being greater 
than unity while they approach unity in a classical limit. Based on the kinematical quanti- 
zation alone, their values are not fixed but subject to quantization ambiguities. Similarly, 
in the gravity sector we have seen a quantum correction function a subject to ambiguities. 
A priori, these quantum correction functions are independent. On the other hand, these 
functions change the dispersion relations of gravitational as well as electromagnetic waves, 
and the corresponding changes in propagation velocities may give rise to concerns regarding 
causality. In particular, the propagation of gravitational waves may become super-luminal 
depending on the precise form of correction functions. 

There are, however, further consistency conditions once the dynamics of the quan- 
tum fields is considered. In a canonical formulation of general relativity, the classical 
constraints Cj form a first class Poisson algebra, i.e. {Ci,Cj} = fl < J (A,E)CK whose co- 
efficients fjj{A, E) can in general be structure functions. The first class nature, i.e. the 
fact that the Poisson brackets of constraints vanish on the constraint surface defined by 
Ci = 0, ensures that the transformations generated by the constraints are gauge and are 
tangential to the constraint surface. Quantum correction functions such as oc(E^) change 
the constraints and thus their algebra. Making sure that the corrected constraints remain 
first class, i.e. that there is no anomaly, provides additional consistency conditions beyond 
those following from the kinematical quantization. As we will see, closure of the corrected 
constraint algebra, in particular for the Poisson bracket of H phen [N] := H G hcn [N] + 
with itself, leads to a relation between all the quantum correction functions in the matter 
and gravity sectors. 

Specifically, the classical Hamiltonian constraint satisfies 



where cross terms between matter and gravity contributions drop out because if em [A] 
couples minimally to gravity. On the other, the gravitational Hamiltonian constraint itself 
satisfies 



where, without loss of generality, we assume the gravitational Gauss constraint to be solved. 
The matter term of expression (jUJ) of Hem [A] 



{HlNx), H[N 2 )} = {HgIN^HgINz]} + {H^lN,}, H EM [N 2 }} 



(61) 



{Hq [Ai] , if g [A 2 ] } = D G [N 1 d a N 2 - N 2 d a Ni\ 



(62) 



{#em[Wi], #em[W 2 ]} = DeuIN^^ - N^N,} 



(63) 



The equations f[6"Tj) . (1621 and fl6"3"|) together thus lead to 



{H[Ni\, H[N 2 }} = D[Nid a N 2 - A 2 <9 a A 1 ] , 



(64) 



where Z}[A a ] is the total diffeomorphism constraint. 
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With quantum corrections we have the gravitational Hamiltonian constraint 

Ht n im = -^gI d * xNa{Et) ^0W\ ^ kFld ~ 2(1 + l2)K l*i) ^ 

which now satisfies 

{Bf'Wi], H% hen [N 2 ]} = D G [a\N^N 2 - N 2 d a N,)] (66) 

(for details see [13]). For the corrected Maxwell Hamiltonian (|55p . on the other hand, we 
have 

{H^tm, HST[N 2 ]} = £> E M[a E M/3 EM (iVia o iV 2 - N^N,)] . (67) 
This can be combined to a first class algebra of the total constraints if and only if 

a 2 = «em/3em , (68) 

such that 

{if phen [iVi], H phcn [N 2 }} = D[a 2 (N 1 d a N 2 - N^Nt)} . (69) 

For linear waves, it is sufficient to use the relation 0681) between the homogeneous parts 
of quantum correction functions, i.e. a 2 = cSem/^em- They appear in the group velocities 



V gw = jJj0A = a and ^EM = y aEM/^EM (70) 

for gravitational and electromagnetic waves. Thus, the requirement of a closed constraint 
algebra, implying (168]) . ensures that there is no violation of causality: the corrected speed 
of gravitational waves agrees with the physical speed of light, which itself is subject to 
corrections. 



7 Wave propagation and lattice refinements of quan- 
tum gravity 

As a further fundamental application, we analyze holonomy corrections in more detail be- 
cause they can give insights into the precise form in which an underlying discrete quantum 
gravity state is being refined during its evolution. Holonomies as multiplication operators 
in loop quantum gravity can create new edges and vertices of a lattice state, and thus 
can dynamically imply its refinements. This can also be described at the effective level 
where, however, the complicated relation to the full theory requires one to refer to several 
parameters describing this refining behavior and in particular the functional form of p,(p) 
used before. Here we show that tensor mode dynamics can be used to restrict the possible 
choices. 
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We parametrize the Hamiltonian constraint as 



TT-phen 
n G 



[N] 



16nG 



d 3 xN 



/sin^r 2 
V ^7 



/ sin /ryfc N 



2^3/2 y ^ 



J k) 



• (71) 



where one parameter is m, the other appears in the power law form jl(p) oc \p\ n . Here 
we have already required that the effective Hamiltonian (1711) has a homogeneous limit 
in agreement with what has been used in isotropic models. This fixes the parameters 
analogous to m in the first two terms to equal one. The parameter for the last term cannot 
be fixed by taking the homogeneous limit and is thus kept free for now. 
Expression (ITTj) provides corrected second order equations 



V J 



(72) 



where 



T< 



Q 



1 

2 

- 2 



sin/17/A", - -x / sin ma jk sin2u7/cY 
(cosm/ry/c — cos Zfi'yk) — I = = I 



M7 
p dp 
p dp 



m/i7 



2/Z7 



v , sin £7* 



smp / yk\ / _ -sinm/ryfc _ -sin/ry/c 

cos fi'jk = cosmp^k 



(73) 



A*7 /V m/i7 /ry 

As before, corrections to the dispersion relation take the form of an effective mass term, 



to ~ 



- I +n>t 
a 



(74) 



where 



Tc 



m„ 



-In 



m 



m 



4 



m 



6 



/x 2 7 2 p \p c 



(75) 



As one can see, this effective mass squared is not guaranteed to be positive for all param- 
eter values. Thus, stability of the perturbation can be used as a criterion to restrict the 
ambiguities. 

One can use anomaly cancellation to relate the free parameters, for which we have 
to refer to vector modes since the tensor mode equations are automatically anomaly-free. 
Specifically, we use the Poisson bracket between the diffeomorphism and Hamiltonian con- 
straints and ensure that it is again linear in the constraints. For simplicity we will consider 
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here only effects of source-free vector perturbations, and correspondingly assume that mat- 
ter constraints vanish. The perturbed Hamiltonian constraint including only vector mode 
perturbations is 



HV hcn [N] 



16nG Js 



d s xN 



V 07 



1 / sin jjrfk N 



2p 3 / 2 \ ^7 



m5E d k 5%) 



+ 



, (76) 



and the perturbed diffeomorphism constraint is 



D G [N a ] 



d 3 x5N c [-p(d k 5K*) - k5 k c {d d 5E%)] . 



(77) 



With the Hamiltonian constraint (1761). we then have 



m/i7 
1 



8ttG 



d 3 xp(<9 c cW J ).4' 



where 



AT 

Vp 



p 



d_ 

dp 



sin iM~ik N 



+ 



sin /ry/c 
07 



k + k 



msin 2/i7&; 



m/ry 



2 sinm/i7A;^ 



(75 



P 



(79) 

The Poisson bracket has terms which cannot be expressed through the constraints unless 
one imposes restrictions on the parameters. To evaluate this, we have to recall that even 
though we write sines in the expression (1711) of quantum corrections, it is to be under- 
stood as a convenient notation to consider leading order quantum corrections. Anomaly 
cancellation up to order k A then leads to the condition 



m 



5 + 2n 



such that 



Tr 



m„ 



(22n 2 - 35n - 5 
V 18 



8ttG 



G*Vp) p 2 



JO) 



11) 



depends on only one remaining parameter n. We have also used the background Hamilto- 
nian constraint to express Tq in terms of the background energy density p. 

The requirement of a positive 'effective mass' squared now implies —0.1319 > n > —5/2, 
restricting the possible functional form of Jx as a function of p. As one can see, some part 
of the otherwise allowed range — l/2<n<0is ruled out here, including a non- refining 
dynamics n = 0. The other limiting case, n = —1/2 of [28], on the other hand, is allowed. 
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8 Discussions 



We have considered tensor mode perturbation equations in Hamiltonian cosmology based 
on Ashtekar variables. In particular, we have derived possible effects of quantum gravity on 
the dispersion relation of gravitational wave propagation in a flat cosmological background. 
Included were typical corrections that one expects from loop quantum gravity, arising for 
inverse volume terms in the Hamiltonian constraint and from the use of holonomies. All 
final results are independent of gauge or other choices in the derivation. 

This shows that inhomogeneities can be considered consistently within a perturbative 
framework of loop quantum gravity. So far, no complete effective Hamiltonian has been 
derived, but several separate effects are known and have at least partially been computed. 
Different types of quantum corrections can thus be studied separately to elucidate possible 
consequences, always keeping in mind that eventually all of them have to be combined for 
a complete picture. The two types of corrections considered here result in rather different 
correction terms in dispersion relations for gravitational waves, which indicates that it is 
reasonable to keep these corrections separate. Typically, only one of them will be dominant 
in a given cosmological regime, and the consequences have different physical consequences. 

Since the magnitude of all the corrections depends on the precise form of a quantum 
state, such properties must be known for a precise quantitative estimate. Qualitative 
implications are, however, clear based on more general principles of loop quantum gravity. 
Also the rate of change of correction terms during cosmic evolution depends on the precise 
state and in particular its refinement. From the tensor mode analysis we have provided 
further evidence that discrete graph states of loop quantum gravity must be refined during 
evolution, supporting the results of [281 EHl HU SSI SB]- Details will also determine the 
precise rate of varying speeds of light and gravitational waves. 

The results provide a viability test of loop quantum gravity already in the absence of 
observations: no violations of causality occur even if quantum corrections in the dispersion 
relations are considered. Along similar lines one has to evaluate more general implications 
of Lorentz symmetries, especially in the context of potential Lorentz violating effects where 
anomaly issues have not yet been considered in the literature. While anomaly calculations 
are difficult for full quantum operators, we have illustrated that partial information can 
be gained economically at the effective level. As seen here, the requirement of anomaly- 
free equations, while allowing for non-trivial quantum corrections, eliminates one effect 
which would otherwise blatantly violate Lorentz invariance. This requires a close relation 
between quantizations of gravitational and matter (especially Maxwell) contributions to the 
Hamiltonian constraint, which is realized by the quantization procedures of loop quantum 
gravity [T9J [31] and tightened by the requirement of an anomaly-free constraint algebra. 
There is thus a weak sense of unification of gravity and matter since quantum corrections 
in the respective terms cannot be independent of each other. 
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